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In the vortex state of a d-wave superconductor, niassless Dirac quasiparticles are scattered from magnetic vortices 
via a combination of two basic mechanisms: effective potential scattering due to the superflow swirling about the 
vortices and Aharonov-Bohm scattering due to the Berry phase acquired by a quasiparticle upon circling a vortex. 
In this paper, we study the superflow contribution by calculating the differential cross section for a quasiparticle 
scattering from the effective non-central potential of a single vortex. We solve the massless Dirac equation in polar 
coordinates and obtain the cross section via a partial wave analysis. We also present a more transparent Born-limit 
calculation and in this approximation we provide an analytic expression for the differential cross section. The Berry 
phase contribution to the quasiparticle scattering is considered in a separate paper. 



I. INTRODUCTION 



The low energy excitations of a d-wave superconductor are the Bogoliubov quasiparticles that reside in the vicinity 
of the four gap nodes at which the order parameter vanishes. The physics of these quasiparticles can be probed 
experimentally via a variety of low temperature electrical and thermal transport measurements. Since quasiparticles 
are part electron and part hole, their energy is well defined but their charge is not. Thus, it is thermal current 
that follows quasiparticle current. Since T ^ Tc, we are well within the superconducting state. Since T ^ Ao (the 
gap maximum), transport is dominated by quasiparticles excited in the vicinity of the gap nodes. (Quasiparticle 
dispersion is therefore given by the anisotropic Dirac spectrum, E = {vjkj + v|fcf )i/2,where Vf is the Fermi velocity, 
V2 is the slope of the gap, and ki and k2 are defined locally about each node. We shall choose our axes such that 
gap nodes are located at zLppx and ^pfY in momentum space.) Furthermore, the temperatures of interest are low 
enough that sources of inelastic scattering are frozen out. This is what we mean by low temperature quasiparticle 
transport. In the mixed state, the remaining energy scales are the impurity scattering rate, Fq, the vortex scattering 
rate, F„, and the temperature. We will focus on the weak-field high-temperature regime where Ty T and T is the 
dominant energy scale. In this regime, the quasiparticles responsible for transport are thermally generated rather than 
impurity-inducecP^ or magnetic field-induced^®'. (Note that thermal transport in the opposite, low T, regime has been 
discussed frequently in the literature^ "^^.) In this high T regime, the physical situation is relatively simple. Thermally 
excited quasiparticles carry the heat current. To understand the thermal transport, we need only understand how 
they scatter from magnetic vortices. Furthermore, since H <C Hc2, the vortices are dilute, separated by distances 
large compared to the quasiparticle de Broglie wavelength. We can therefore learn a lot by considering the scattering 
of quasiparticles from a vortex and assuming that all such scattering events are independent. We take the following 
path. In Sec.|IT] we develop our picture of a nodal quasiparticle scattering from a vortex. Then in Sec. |III[ we consider 
the Bogoliubov-de Gennes equation, make a singular gauge transformation, and shift the origin of momentum space 
to the location of one of the gap nodes. The resulting problem is one of an (anisotropic) Dirac fermion scattering 
from an effective non-central potential (due to the superflow) in the presence of antiperiodic boundary conditions 
(a consequence of our gauge choice). After obtaining the quasiparticle current functional in Sec. IV we define our 
vortex model in Sec. [V] and make several important approximations. As far as scattering is concerned, the vortex has 
two parts: a circulating superfiow and a Berry phase factor of (-1) acquired upon circling a vortex. In this paper, 
we consider scattering due to superfiow current in a single vortex, without the Berry phase effect. (We consider the 
Berry phase contribution in a separate papei^.) Furthermore, we neglect the anisotropy of the Dirac dispersion and 
take Vf = V2- Although these approximations will prevent our results from being quantitatively accurate, we argue 
that our qualitative results refiect the essential physics of the problem. This notion is supported by Ref. [131 wherein 
this model, with these approximations, was used to derive a thermal Hall conductivity in good agreement with the 
qualitative features of the experimental results for YBa2Cu30e,gg by Zhang et a^^. In Sees. VI and VII we expand 



the quasiparticle wave function in angular momentum eigenstates and obtain the resulting radial equations. We solve 
these equations outside the vortex in Sec. |VIII| and inside the vortex in Sec. |IX[ In Sec. |X] we match solutions and 
sum over all four gap nodes in order to obtain the single vortex scattering cross section. Results of our numerical 
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Figure 1: Schematic depiction of vortices in the mixed state of an extreme type II superconductor. Vortices are pinned to 
randomly distributed pinning centers and are separated by an average inter-vortex distance, 2R. Since ^ ^ R <^ X, vortex cores 
are point-like while the magnetic field profile of each vortex is quite extended. Therefore, the magnetic field is approximately 
uniform across the sample even though the vortices are well separated. 



calculations are presented in Sec. XI In the same section (Sec. XI I we write a closed-form expression for the differential 
cross section using the more transparent Born approximation details of which are given in Appendix ^ Conclusions 
are discussed in Sec XII Note that the results of this calculation were used in Ref . il3j (without the detailed derivation 
presented here) to calculate thermal transport coefficients. 



II. SINGLE VORTEX SCATTERING OF BOGOLIUBOV QUASIPARTICLES 



Consider the mixed state of a d-wave superconductor. Our picture is that of Fig. [T] In the presence of a magnetic field 
{H > Hci), magnetic vortices penetrate the sample (a 2D Cu02 layer). Vortices are distributed randomly, pinned to 
local defects. The cuprates are extreme type II superconductors in which the coherence length, ^, is much smaller 
than the penetration depth, A. As a result, while the vortex cores may be well separated, the magnetic field profiles 
overlap significantly such that there is little variation in the magnetic field across the sample. We therefore adopt the 
extreme type II limit of ^ and A — ^ oo and take the magnetic field to be constant, H = Hz. In this limit, there 
are only two remaining length scales. The first, 1/k, is set by the temperature such that k = E/hvf = ksT/fiVf. For 
anisotropic Dirac nodes, k — \/kl + ^ll^s/f/)^. Note that in the isotropic limit, this reduces to the magnitude of the 
quasiparticle momentum (measured from a node). The second length scale, R, is half of the average distance between 
vortices. With one flux quantum per vortex, HttR? = $o = hc/2e, so we define R = yJfic/eH. In terms of R, we can 
define the (2D) density of vortices to be Uy = H/^q = 1/ttR^. Note that the ratio of these two length scales yields 

kR^'l^.[^=jL (1) 
which is the inverse of the argument of the Simon-Lee scaling functionJi^], 

We wish to consider a quasiparticle, in a state with a particular energy and current, scattering from a magnetic vortex 
into another state, of the same energy, with a different current. At low temperatures, quasiparticles are excited in the 
vicinity of the four gap nodes in momentum space. The nodal structure of the Brillouin zone is depicted in Figj2j^b). 
Following Simon and Lee^**, we have defined coordinate axes, rotated by 45° from the usual crystal axes, such that the 
gap nodes are located on-axis at ±pj?x and ±ppy. In the neighborhood of a node, the quasiparticle excitation energy 



is given by the anisotropic Dirac spectrum, Ek — y ^/^i + "^2^2- where the fci and k2 axes are defined, respectively, to 
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Figure 2: Schematic depictions of momentum space and coordinate space, (a) Quasiparticle scattering from a single vortex. 
An incident plane wave scatters from an origin-centered vortex into an outgoing radial wave.(b) Gap nodes in tlie Brillouin 
zone of a Cu02 layer. Ellipses denote surfaces of equal energy. Arrows represent the quasiparticle group velocity in the vicinity 
of each node. 

be perpendicular and parallel to the local Fermi surface. The surfaces of equal energy are therefore ellipses centered 
about each of the nodes. The quasiparticle current is given by the group velocity = Vfc-Bfe, directed outward from 
the node centers, perpendicular to the ellipses of constant energy. 

In coordinate space, the quasiparticle is represented as an incident plane wave. The single (/ic/2e) vortex is centered 
at the origin and surrounded by a circulating supercurrent. Due to the interaction between the quasiparticle and the 
vortex, the quasiparticle is scattered into an outgoing radial wave. The situation is depicted schematically in Fig.|2ja). 
By solving the Bogoliubov-de Gennes equation for the quasiparticle wave functions and considering the nature of the 
quasiparticle current before and after such a scattering event, we will compute the quasiparticle cross section for 
scattering from a vortex in a d-wave superconductor. Our analysis shall be similar in spirit to that conducted for 
the s-wave case by Clear jff^EH^ -^j^q extended the work of Caroli, de Gennes, and MatricoiJ^^ (on vortex core bound 
states) to the problem of quasiparticle scattering from a vortex in an s-wave superconductor. 



III. BOGOLIUBOV-DE GENNES EQUATION 

Consider the Bogoliubov-de Gennes (BdG) equation for a d-wave superconductor in the presence of a constant per- 
pendicular magnetic field, A — ^Hrcj), and with an order parameter that winds once about the origin, A(r) = Aoe*"^: 



H' = 



A'* 



A' 



(2) 



1 

2m 



P--A -E_ 



(3) 



A' = ^{p;,{p;,A(r)}} 
Pf 



Ap 



2 A(r)(a,9,0) 



(4) 



Here p — —ihV, {a, 5} = {ab + ba)/2, E is the quasiparticle energy, and the form of the gap operator is that required 
to provide d-wave symmetrj^l^IIll Entangled within the complex differential form of the gap operator. A, is the fact 
that upon circling an hc/2e vortex, the quasiparticle acquires a Berry phase factor of (-1). The Hamiltonian would be 
simplified if we could effectively strip the gap function, A(r), of its phase. This is accomplished via the application 
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of the singular gauge transformation 



U = [ ""J^^^' J/, ] $(r) = U-'^{r) H = U-^H'U. (5) 



In this gauge, known as the Anderson gauge, 



i7$ = E<^ (6) 



^ = ^3 [(P + T-aP.)' - pI] + n [2p^py\ (7) 

2pF 2pF 



where 



P.M. »v*-Sa. 5(1-1,)^ ,s, 

is the gauge invariant superfluid momentum (superflow), Vf = pf/tti, V2 — Aq/pf, and R = y/hc/eH . Although 
this form of the Bogoliubov-de Gennes equation is much simplified, the boundary conditions have become more 
complicated. While the original wave function was defined with periodic boundary conditions, ^'(r, (jj) — ^(r, + 27r), 
the transformed wave function is not single-valued and has antiperiodic boundary conditions, <i>(r, (j)) = — $(r, (j) + 2Tr). 
Hence, we have introduced a branch cut such that with each trip around the origin, the wave function changes sign. 
In effect, the Berry phase contribution has been extracted from the Hamiltonian and encoded in the antiperiodic 
boundary conditions imposed on the wave function. Note that when considering a system of many vortices, the 
resulting sea of branch cuts can be problematic. In re sponse, Franz and Tesanovic introduced an alternate gauge 
transformation better suited to the vortex lattic^^^ESI. However, for the single vortex case that we consider, the 
Anderson gauge described above is most convenient. All observables, such as differential cross section or transport 
coefficients, are independent of gauge choice. We would like to emphasize that the Berry phase effect is not an artifact 
of our gauge choice. It is physical, and an intrinsic part of the original Bogoliubov de Gennes equations. In the 
untransformed Hamiltonian, which has periodic boundary conditions, it resides in the complex differential form of the 
gap operator. The gauge transformation drastically simplifies the Hamiltonian at the cost of introducing antiperiodic 
boundary conditions. Thus, effectively, the Berry phase contribution is extracted from the Hamiltonian and moved 
to the boundary conditions. 

Recall that the low energy excitations of a d-wave superconductor are concentrated about the four nodal points in 
momentum space where the gap vanishes. We can therefore further simplify our Hamiltonian by shifting the origin 
of momentum space to the location of one of the nodes. Shifting to node 1 



we find that 



Px ^ PF + Px Py ->■ Py (9) 

H = Hd+Hc (10) 

Hd = Vf [pa:T3 + apyTi + Psx] (11) 

Hc^^ [{p" + Ps)rz + 2P, • p + a2p.,pyTi] (12) 

IpF 

where a = V2/vf and we have used the fact that Ps — Ps{r)(j} to commute p with P^. Here H is written as the sum 
of a linear (Dirac) Hamiltonian, Ho, and a quadratic (curvature) Hamiltonian, He, to emphasize that the second 
term is smaller than the first by a factor of E/Ef- We will focus on the dominant term, Hf, , in the remainder of 
this paper. In what follows we shall seek solutions to the Bogoliubov-de Gennes equation developed above. However, 
we must first consider the nature of quasiparticle current in a c?-wave superconductor. 
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IV. QUASIPARTICLE CURRENT 



In order to calculate the cross section for a quasiparticle scattering from a vortex, we must be able to write down the 
currents corresponding to the incident and scattered wave functions. Since the incident and scattered currents will 
be considered in the far field where the quasiparticles are free, we wish to determine the quasiparticle current as a 
functional of $ for P, = 0. Setting P, = in Eq. ^ we find that the BdG Hamiltonian becomes 



H 



He A 

A* -h: 



He 



2p 



V2_ 
'2pF 



2dj:dy 



(13) 



Following Ref. '501, we can write down a mean-field Hamiltonian in terms of the real-space electron cre- 
ation/annihilations operators, and '^a- 



Heff = / d 



t,T,t 



(14) 



Since quasiparticles carry well-defined spin (and heat) the quasiparticle current is equal to the spin current divided 
by the quasiparticle spin. The spin density operator is expressed as 



where s = 1/2 and the spin current operator is, in turn, determined via 

V.f(r) = V=z[p^(r),i7,//]. 
Making use of the fermionic anticommutation relations, 

{*„(r), ^p{r')} = {vl/t^(r), *t (r')} = {^l{r), ^p{r')} = <5„^<5(r 



we find that 



(15) 



(16) 



(17) 



(18) 



where rja = ±1. Noting that each of these terms can be manipulated into the form of a divergence and making use 
of Eq. (16 1 we can write the spin current operator as 



2ipF 



(19) 



where V = + dyj and Vj^ = 9j,x 4- d^y- The real-space creation/annihilation operators, vj/t and ^^^i can be 
expressed as a weighted sum of Bogoliubov operators, and 7„q, via 



(20) 



where it„ and Vn are particle and hole wave functions for state n. Plugging these forms into the spin current operator 
and evaluating for state i and spin a yields an expression for the spin current as a functional of and vg. 



{£a\-f\icT) 



srja 
2ipF 



[Vfiul'Vug — W^Vtif) -I- f2(w|V_LW£ + f|V_LU£) — c.c. 



(21) 



Dropping the eigenstate label, dividing by the spin sr]^, and writing the result in terms of the particle-hole 2- vector 



$ = 



$t ^ {u*,V*) 



(22) 
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we obtain the quasiparticle current 



j = — Ini [w/^VaV* + V2^^ TiV _i<S>] . (23) 
Pf 



As in the preceding section, it is convenient to shift the origin of momentum space to a nodal point. Shifting to node 
1 yields 

3=3D+jc (24) 

^Vf^^Tsk + aTiy)^ (25) 



^/ T 

Jc = — Im 

PF 



(26) 



where a — V2/vf. Once again, the second (curvature) term is smaller than the first (Dirac) term by a factor of E/Ep 
and shall be neglected in what follows. 



V. MODEL AND APPROXIMATIONS FOR SINGLE VORTEX 

Given the BdG Hamiltonian and the quasiparticle current functional, we are ready to consider a model of quasiparticle 
scattering from a single vortex. In reality, each quasiparticle encounters a sea of vortices separated by distances on 
the order of 2i? = yj Ahc/eH. The superflow circulating around each vortex is peaked near the vortex centers but 
falls off slowly and overlaps in the regions between vortices. The total flux through the sample is equal to one (/ic/2e) 
flux quantum per vortex but is distributed smoothly since the penetration depth is much larger than the inter- vortex 
distance. In order to model this situation via single vortex scattering, we approximate the effect of neighboring vortices 
by cutting off the superflow distribution about our single vortex at a distance R from its center. By construction, the 
flux through this circle is exactly one flux quantum. Thus, the superfluid momentum takes the form of Eq. ([8| inside 
the vortex (r < R) and is zero outside the vortex {r > R) 



V,{v) = U\-^]0{R~r)4^. (27) 



Note that at the vortex edge {r — R), the phase gradient and vector potential terms cancel and the superflow goes 
continuously to zero. We shall consider a quasiparticle plane wave incident upon such a vortex and calculate the 
scattering cross section. Before going forward, however, we make the following additional approximations. 



First, recall from Sec. |III| that, with our gauge choice, the wave function must obey antiperiodic boundary conditions 
such that it changes sign each time the quasiparticle winds around the vortex. This is the Berry phase contribution to 
the problem and it has the interference effect of an Aharonov-Bohm tt-Aux at the origin. Thus, even in the absence of 
a superflow, it is possible for quasiparticles to scatter, due solely to the antiperiodic boundary conditions. The effect 
of this Berry phase is considered in a companion paper. 

Second, recall that the quasiparticle excitation spectrum for a d-wave superconductor is that of an anisotropic Dirac 
cone with two characteristic velocities, Vf and V2, which can be quite different. This anisotropy enters the BdG 
equation and current functional in the form of the parameter a — V2/vf. The fact that a is not equal to one has the 
effect of complicating the form of the BdG equation in polar coordinates. Thus, in order to make this partial differential 
equation separable, we should scale out the anisotropy. This can be accomplished by scaling the y-coordinate by a 
factor of a. In the scaled coordinates, the elliptical gap nodes become circular and we obtain the isotropic Dirac 
equation. Unfortunately, this scaling also has a negative consequence. Whereas the vortex was originally circular in 
real space, the scaling makes it elliptical. Thus, in the scaled coordinates, the superfluid momentum (and the vortex 
boundary) become much more complicated. While such complications can be dealt with (via a significantly more 
involved computation), they were not considered in this investigation. Possible consequences of anisotropy have been 
studied in Ref. ? , where it was noted that there are no drastic experimental signatures due to the anisotropy of the 
Dirac cone. Thus, for simplicity, rather than scaling out the anisotropy, we consider the more straightforward case 
where the gap nodes are isotropic and Vf = V2- For this isotropic case, both the gap nodes and the vortex are circular 
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and we can separate the BdG equation in polar coordinates. However, since Vf exceeds V2 by a factor of 10 to 20 in 
the cuprates^iJ this is clearly an approximation. As a result, we expect only qualitative agreement with reality. That 
this is the case was demonstrated by Durst, Vishwanath, and Lee^^who used this approximation to provide a clear, 
though qualitative, explanation of the thermal Hall conductivity measurements of Ong and co-workerd^ 

Finally, recall that the BdG Hamiltonian can be expressed as the sum of a linearized (Dirac) part. Ho, and a quadratic 
(curvature) part. He- We see that He is small compared to Hjj as long as we are sufficiently far from the vortex center 
(r > 1/pp). Hence we must cutoff our model at the scale of the vortex core ^ 10/pp). As our final assumption, 
we select a reasonable core size and model the vortex core as a region with vanishing superflow. We now have a 
well-defined scattering problem, which is solved considering the linearized Hamiltonian. Since we shall only consider 
the isotropic case from this point forward, the Hamiltonian becomes 

Hd = Vf [tsPx + Tipy + Psx] (28) 

the quasiparticle energy is 

E^Vf^kl + kl=Vfk (29) 

and the quasiparticle current functional takes the form 

ji5 = w/$^(T3A + Tiy)$ (30) 

Note that the effective potential induced by the superflow is smooth on the scale of 1/pp. The Fourier components for 
large-momentum inter-node scattering are therefore small. Hence, we expect that the contribution to the scattering 
cross section of inter-node scattering will be subdominant to that of intra-node scattering. Thus, for simplicity, we 
include only intra-node scattering in what follows. 



VI. ANGULAR MOMENTUM EIGENSTATES 



Consider the linearized Hamiltonian, Hpi, in the isotropic limit, Eq. (28). Note that for Pg = 0, this is just the Dirac 
Hamiltonian for massless spin- 1/2 fermions in two dimensions 

H%^Vf[T3p,+TiPy]=Vf[a-p + l3m] m-0 /? = 7° a'' = 7V (31) 
with the 7-matrix representation 

7° — T2 — iTl 7^ = — «T3 — > (3 = T2 — T3 O? = T\ (32) 

For the physical Dirac equation we know that the total angular momentum operator takes the form, J = i + |l] 
where E — \tj,jkl^^^ ■ Therefore, in the above 7-matrix representation 

Evaluating the commutator with H'j^ we find (as expected) 

[J, Hl\ = -T2K + = 0. (34) 

Therefore, there exists a complete set of simultaneous eigenstates of i/^ and J. 
The eigenstates oi L — —i-^ take the form 

L* = £* ^ VP ^ e^^^ ) (35) 
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and the eigenstates of E = r2 take the form 

E* = A* ^- A = ±l 



c(r, (/)) 



Therefore, the simultaneous eigenstates of J and have the form 
J$„ = (n + l/2)$„ ^ $„ = /„(r)e"'^ 



9n{r)e 



1 

— Z 



1 

— i 



(36) 



(37) 



where the radial functions, /„(r) and gnif), are determined from the solution of a pair of coupled radial equations. 
Since we have neglected the Berry phase contribution and adopted periodic boundary conditions for this case, the 
requirement of single- valued wave functions demands that n is an integer. 

For Ps 0, Hd looks like the massless Dirac equation in the presence of an effective scalar potential, V — VfPgx — 
— f/Ps(r) sin0. Since the effective potential is non-central, it mixes angular momentum eigenstates. (See Refs. 
and EH for a discussion of electron scattering from a similar (dipole) potential.) While the general solution can still 
be expressed as a linear combination of angular momentum eigenstates 



9nir)e 



»(n+l)0 



1 



(38) 



the radial equations for different n are now all coupled together. Nonetheless, general solutions can still be written in 
the form of Eq. (38 1. In the following section, we apply Hu to a wave function of this form and proceed to determine 
the resulting radial equations. 



VII. RADIAL EQUATIONS 



We shall now plug the general form of our wave function, Eq. (38 1, into the Bogoliubov-de Gennes equation. For this 
purpose, it is useful to write the BdG Hamiltonian, H = Hd, in polar coordinates. Doing so, we find 



H 



D = -IVf 



. „, 9 (t ■ f ) . d „ . . 

(r • r) — H «T2^ - Psismc 

or r o4> 



(39) 



where (t • f ) = T3 cos0 -1- ti sin0 and we note that (t • f)^ — 1. It is also useful to write the wave function as 

a>(r) = ^[/„(r)x++5n(r)x,T] 



(40) 



where we have defined 



„i(n+l)0 



1 



(41) 



Note that (r • f ) transforms xt, i^^o Xn y'lce versa. 



{t ■ f )Xrt ^Xn {t- T)Xn = xt 



(42) 



and also that cosines and sines shift the x's up and down in angular momentum. 



, ± Xn+l "t" Xn-1 . I ± Xn+l Xn-1 
COS</)X„= ^ S\TL(t)Xn = 77-- ■ 

2 2,1 



(43) 



Making use of these relations, we see that both and i?$ can be written as sums over Xn weighted by coefficients 
that are functions only of the radial coordinate. Equating coefficients of Xn both sides of the BdG equation and 
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dividing out an overall factor of —iv/k yields a set of fully coupled differential equations for /„ and 

P. 



dfn 

dp 



In 



p 



n-1 — ffn+lj 







(44) 



dgn 
dp 



p 2 



= 



(45) 



where we have defined a dimensionless radial coordinate, p = kr, and a dimensionless superfluid momentum, Ps{p) = 
(l/p — p/ {kRY) /2. In this manner, the solution of the BdG equation is reduced to the that of a system of coupled 
ordinary differential equations. 



VIII. OUTSIDE THE VORTEX 



According to our model of the single vortex, the superfluid momentum vanishes at a distance, R, from the origin. Thus, 
for p > fci?, we sit outside the vortex. Here Pg = and the quasiparticles are free. Since the superfluid momentum 
(Eq. [8| includes the effects of both the order parameter phase gradient and the magnetic field, quasiparticles outside 
the vortex are subject to neither. Thus, for p > kR, we consider quasiparticles subject only to the free Dirac 
Hamiltonian, H^, and with current defined via the linearized current functional, In this regime, it is possible 

to obtain explicit solutions to the BdG equation, define incident and scattered wave functions, and, by constructing 
these incident and scattered waves from the free basis functions, write down an expression for the scattering cross 
section in terms of the coefficients of the basis functions. Then we need only match solutions with those inside the 
vortex to obtain the cross section. 



A. Free Solutions 



For Pg = , the radial equations, (44 1 and (45 1, take the dramatically simpler form 

d n + 1 ' 



d n 
dp p 



In ^9n 



dp 



gn = ifn- 



(46) 



Note in particular that, while /„ is coupled to g„, the equations for functions of different n are independent. Elimi- 
nating gn from the equations above yields 



1 



dp"^ p dp \ p 
which is the defining equation for the Bessel functions. Therefore, 

fn{p) = AnJnip) + BnY„{p) 



(47) 



(48) 



where A„ and B„ are complex constants and Jn{p) and Yn{p) are Bessel functions of the first and second kind. 
Substituting back for gn then yields 



gnip) = iAnJn+l{p) + i-B„F„+i(p) 



(49) 



where we have used the Bessel function identity {d/dp — v/ p)Zi, = — Hence, we can write down the free wave 
function 



$(r) = ^ {AnJn + BnYn)e"^^ (])+ i(A„ J„+i + B„y„+i )e'("+i'^ ( ^. 
where the coefficients shall remain undetermined until we match with solutions inside the vortex. 



(50) 
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B. Incident and Scattered Waves 



We wish to obtain the cross section for a plane wave, with quasiparticle current in the incident direction, scattering off 
a vortex as a radial wave, with quasiparticle current in the scattered direction. If the incident momentum is k = (fc, 9) 
and the final momentum is k' = (fc, (p), then the incident direction is the direction of the group velocity at momentum 
k and the scattered direction is the direction of the group velocity at momentum k'. For general, anisotropic nodes, 
the group velocity need not be parallel to the momentum. However, for the isotropic case that we consider 

VG(k) = = Vf^x + V2^y = Vf (cos 6lx + sin 6ly) = Vfk (51) 

OK I^k Ek 

and the group velocity and momentum are parallel. Therefore, if denotes the incident wave function and denotes 
the scattered wave function, then we require 

(cos6lx + sin6'y) - k - (cos ^x + sin (/)y) - k' - f . (52) 



Recall that outside the vortex, quasiparticles are subject to neither an order parameter phase gradient nor a magnetic 
field. Thus, the incident wave function is a plane wave. This is consistent with the well-known results of Franz 
and Tesanovic (Franz and Tesanovic, PRL 84, 554 (2000)) who showed that the low-energy quasiparticle states of a 
d-wave superconductor in the vortex state are Bloch waves of massless Dirac fermions rather than Landau Levels. 
(For a discussion of the analyses that led to this important result, the reader is referred to Refs. [25J|ini[2Zl andlT^V 
Inspection of the form of the current functional, jjj = Vf^^T^x + Tiy)$ reveals that the appropriate incident plane 
wave is 



(53) 



smce 



sm - X 

2. 



2 sin - cos - I y 
2 2 ' 



Vfk. 



(54) 



Note also that this form solves the BdG equation, as it must in the absence of the vortex. The appropriate scattered 
radial wave is then 



since 




Here f{(f> — 6) is the scattering amplitude and the e*"^/^ prefactor has been added to make the wave function single- 
valued. 



C. Constructing the Cross Section 



We shall now construct, from our free solution basis functions, an asymptotic wave function containing the correct 
incident and scattered waves. The factor in front of the scattered wave will then be our scattering amplitude. 



We begin with the free wave function obtained in Eq. (50 1. Note that in the asymptotic limit, the integer-index Bessel 
functions can be expressed as 



Jn{p) 




cos(p-^/4- |n|7r/2) Y,,{p) 




sin(/9 - 7r/4 - |n|7r/2) 



(57) 
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where ?]„ = 1 for n > and rjn = (—1)" for n < 0. Therefore, if we shift n ^ n — 1 in the second term of Eq. (50 1 
plug in the asymptotic forms, and decompose the sines and cosines into exponentials, we obtain 



TTp 2 



(58) 

We can replace our two complex constants. An and i3„, with two new complex constants, a„ and 6„, by defining 

An - iBn = z"e-^("+5)'' (1/2 + e^l 6„) (59) 



An + iB„ = z"e-^("+5)'' (1/2 + (-l)"a„) . 
Making use of these definitions, noting that 77„i" — i'"', and reorganizing terms, we find 



2 1 

cos(p — 7r/4 — |n|7r/2)- 
np 2 



2 1 
7rp2 



Noting the Bessel function expansion of a plane wave 

g»pcos(0-e) ^^i«j„(p)e™(*-'') =^i"e"(^-^)?7„yjcos(p- ^4- |n|7r/2) 



and shifting n — > n + 1 in the final term of Eq. (61 1 yields 



sm 2 



I cos ■ 



Regrouping terms and defining tp = (t> — 0, this becomes 



$ = e 



:k-r 



COS : 



e'2 



6^'='^ / COS ^ 



where 



(60) 



(61) 



(62) 



(63) 



(64) 



(65) 



5(<P) 



(66) 



Since we have yet to restrict a„ and 6„, this asymptotic wave function is still totally general. However, we have 
succeeded in rearranging it into a suggestive form. The three terms above are easily understood. The first is our 
incident plane wave, the second is the outgoing radial wave, and the third is an incoming radial wave. By construction, 
we require an incident plane wave and an outgoing radial wave. To realize this scenario, our asymptotic boundary 
conditions require that there be no additional incoming wave. Thus, we require g{ip) = and must therefore set 
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a„ = for all n. In terms of our original constants, this restriction requires that 
and sets 



(67) 



b = e-^J f — - 1 



Furthermore, we obtain a simple form for the asymptotic wave function 



$(r) = e 



sm 7j I Vr V sm : 



(69) 



where we recall that k = (fc, 9), r = (r, 0), and ip = (p — 6. 

Applying the current functional, jo, to both the incident and scattered parts of this wave function yields the incident 
and scattered current density 



Then the differential cross section is 



= I/Ml' 



Integrating over (p yields the total cross section 

By weighting with a factor of (1 — coscp), we obtain the transport cross section 

(7|| = J dip ^(1 - cosy?) = ^ y^b„ [b„ - {bn+i + 6n-i)/2] 

and by weighting with a factor of sin (p we find the the skew cross section 

aj_= [ dip simp = ^y^^n [bn+i - &n-l]* • 



(70) 



(71) 



(72) 



(73) 



(74) 



All that remains is to calculate the coefficients 6„. These, of course, are determined by the details of the quasiparticle 
scattering. Hence, we must now look inside the vortex. 



IX. INSIDE THE VORTEX 



Inside the vortex, the situation is more complicated. For p < kR, the superfluid momentum is nonzero and takes the 
form, Pg = (1/ p— p/ (kR)^) /2. Therefore, even the linearized Hamiltonian, Hjj, mixes angular momentum eigenstates. 
As discussed by Simon and LefP, Y(P, and VishwanatlP, since the linearized Hamiltonian is time-reversal invariant, 
the resulting skew cross section (after summing over all four nodes) is zero. We solve for the radial functions, fn{p) 
and gn{p)- The radial equations look as follows: 



dfn 
dp 



-fn + ifjn + -^{9n-l " ffn+l) 



P 







(75) 
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dp 



9n + «/n + -7rUn-l - Jn+l) 

P 2 



= 



(76) 



Since the radial equations are fully coupled by the superfluid momentum terms, the equations for all n should be solved 
simultaneously via an infinite-dimensional matrix equation. By numerical necessity, we shall cut off the coupling at 



some large n such that we consider a total of N angular momentum eigenstates where 



N 



< n < 



N 



1. This is 



physically reasonable as we do not expect very large angular momenta to have a significant effect on the low energy 
physics. The result is a 2N x 2N matrix differential equation 



dz 

dp 



M(p)z 



(77) 



where z(p) is a 27V-component vector containing the /„'s and the (7„'s, M(p) is a 2N x 2N matrix. Since only 
neighboring angular momenta are coupled, M{p) is a rather sparse matrix. Furthermore, due to the form of Eqs. (751 
and ([76| and the simple p-dependence of the superfluid momentum, this matrix can be written as 



M{p) = B 



1 
P 



Aip 



(78) 



where B, Aq, and Ai are constant 2N x 2N matrices. In terms of this matrix notation, our procedure will be to solve 
the homogenous equation. 



A. Homogeneous Solutions via Method of Probenius 



Consider the homogeneous equation 

dz, 

dp 



M{p)z 



dz 

p— = (B + Aop + Aip2)z 



(79) 



Since M(p) diverges as 1/p at the origin, the equation has a regular singular point at p = 0. Thus, its solution 
requires local analysis about the origin. In particular, we shall employ a matrix generalization of what is known as 
the Method of Frobenius^^. 



Consider a solution of the (Frobenius) form 



a-\-m 



(80) 



m=0 



where a is a complex number and is a vector coefficient for each integer m. Plugging this into our differential 
equation yields 



[(B - (a + m)) a,„p"+" + Aoa™p"+"+i + Aia™p"+"+2] = 0. 



(81) 



Shifting m — m — 1 for the Aq term, shifting m — )■ m — 2 for the Ai term, and noting that the total coefficient of 
each power of p must equal zero, we obtain the following equations: 



(B - a) ao = 



(82) 



(B - (a+ l))ai = -Aoao 



(83) 



(B - (a + rn)) a,,„ = -Aoa,„_i - Aia„_2 m ^ 2, 3, 4, . 



(84) 



Note that Eq. (82) is just the eigenvalue equation for B. It is solved if a is an eigenvalue of B and ao is the 
corresponding eigenvector. This eigensystem has 2N solutions which we can label ttfc and ag for fc = 1 . . . 2N . The 



other two equations, (83 1 and (84l, define matrix recursion relations. If B— (afe+m)! is non-singular for m — 1, 2, 3, 
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then we can solve these equations to obtain given a§ and ak- This condition will always be satisfied for the largest 
eigenvalue of B but will fail for any ak for which ak+mis also an eigenvalue of B. In general, this circumstance would 
require additional solutions not of the Frobenius form. Fortunately, for Pg 7^ 0, the eigenvalues of B are separated by 
non-integer numbers. Therefore, we can obtain 2N solutions of the Frobenius form 

z'' = J2 a^P^'^"' fc = 1 . . . 2iV (85) 

m=0 

and use the equations above to calculate the vector coefficients recursively. Due to the constraints of our numerics, 
we cutoff each series at m = mmax- As such series represent expansions about the origin, the more terms we include 
the larger the p at which the solutions will be valid. Since the differential equation contains no singular points besides 
the origin, the solution can be pushed as far from the origin as necessary by including additional terms. Since our 
solutions need only be valid out to the edge of the vortex, p = kR, we shall cutoff each series at the mmax for which 
the solutions are valid out to kR. 



B. Boundary Condition at Origin 



Although all 2N solutions, z , satisfy the homogeneous equation (79), we must reject solutions that show unphysical 
behavior at the origin. While it is permissible for the wave function or the current density to diverge at the origin, we 
must require that observable quantities like the total probability at the origin and the total current passing through 
the origin be well behaved. Since solutions of the Frobenius form yield contributions to the wave function of the form 
<i>(p — > 0) ^ for small p, we seek conditions on the allowed values of a. 

Consider a circle of radius e — )■ about the origin. The total quasiparticle probability within the circle is 

Prob= f 0)$(p-^ 0) - f {pdp)p"p°'r^ f dpp'^°'+^. (86) 

J Jo Jo 

If the origin is to be just like any other point in space, then this probability must vanish as e — > 0. Hence, we 
require a > — 1. While this restricts the allowed solutions, it turns out that the quasiparticle current provides a more 
restrictive condition. 

Now consider a semicircle of radius e — > about the origin, oriented about the direction. The total current passing 
though such a semicircle in the 6 direction is 

Ie= / (ed0)w/$^(p^O)(T3X + riy)$(p^O)-^-ee"e"~e2a+i^ (87) 
For the origin to be physical, this current must vanish as e — ?> 0. Hence, we require a > —1/2. 

For Pg ^ 0, the eigenvalues of B are always such that N of them are strictly larger than —1/2 and N of them are 
strictly smaller than —1/2. Thus, in order to satisfy these boundary conditions at the origin, we shall keep only the 
z''', for which > —1/2. The full solution to the homogeneous equation is then 

N 



=IIcfcz'=(p) 



k = l 

where the Ck are N complex coefficients to be determined via solution matching at the edge of the vortex. 



X. FOUR-NODE CROSS SECTION 



For p > kR, we found 2N radial solutions (Bessel functions) and reduced the 2N coefficients (yl„,i?„) to N coefficients 
with the condition that there be no incident radial wave in the asymptotic limit. (This required i?„ — i{An — ^°).) 
For p < kR, we found 2N radial solutions (Frobenius series) and reduced the 2N coefficients to N coefficients with 
the condition that the solutions be well-behaved at the origin. (This required = for < — i.) At p = kR, all 2N 
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radial solutions must be continuous. (Note that since the Dirac equation is 1st order, only continuity is required, while 
for the 2nd order Schrodinger equation, differentiability would also be required.) Thus, we have 27V equations which 
can be solved for the remaining 2N coefficients by matching solutions at p = kR. We solve the resulting 2N x 2N 
matrix equation via LU decomposition to obtain Ck for 1 < fc < A'^ and An for — y <n< y — 1- 

As discussed in prior sections, solutions are matched to obtain to obtain the coefhcients to get scattering cross sections. 
Recall that at the outset of this calculation, we shifted the origin of momentum space to the center of node 1. Thus, 
in the discussions that followed, we have been considering quasiparticles scattered from one state in the vicinity of 
node 1 to another state in the vicinity of node 1. The resulting cross section is therefore only the cross section for 
these node-1 quasiparticles. However, given a quasiparticle current in any particular direction, quasiparticles from all 
four nodes will contribute equally. Thus to obtain the physical cross section, we must average over the cross sections 
for quasiparticles at each of the four nodes. 

Our results for node 1 can be easily generalized to node j = {1, 2, 3, 4} by transforming coordinates to those appropriate 
to node j. In accordance with the d-wave structure of the gap, we can define a local coordinate system at each of the 
four nodes with a ki axis pointing along the direction of increasing and a k2 axis pointing along the direction of 
increasing A^. Note that while nodes 1 and 3 define right-handed coordinate systems, nodes 2 and 4 define left-handed 
coordinate systems. We can therefore transform from node 1 to node j simply by rotating our incident and scattered 
angles {6 and (f) and then changing the sign of these angles to account for the handedness of the local coordinate 
system. 



Node 1 
Node 2 
Node 3 
Node 4 



^1 = W (Pi^ (p ipi ^ (pi - Ui ^ if 

% = -{Q-l) 02 = -(0-f) ^1^^2-02 = -^ 

9^ = 9 + n 03 = + TT (^3 03 - (?3 = 1^ 

94 = -(^+f) 04 = -(0+f) (^4-04-^^4 = -'/' 



(89) 



Thus, to obtain results for quasiparticles about node j, we need only input each 9j and take the output as a function 
of {—iy~^^(p. Then the physical cross sections are 



da 

dip 4 



J = l ^ .7=1 



Note that we do not account for quasiparticles that are scattered from one node to another. However, since the effective 
potential (induced by the superflow) is smooth on the scale of 1/pf, the Fourier components for large-momentum 
inter-node scattering will be small. Hence, it is reasonable to neglect such contributions. 



XI. RESULTS 



The procedure outlined in the preceding sections was implemented numerically for a range of intervortex distances, 
kR = /\/H, from 0.5 to 15. The larger the kR, the more angular momentum eigenstates that contribute to the 
vortex scattering cross section, and the more runtime-intensive the computation. For kR = 15, we set iV = 46 and 
considered total angular momenta, n + 1/2, ranging from -22.5 to 22.5. Calculating solutions to the BdG equation 
both inside and outside of the vortex, the program matches solutions at p — kR to obtain the coefficients, b„, first for 



node 1 and then for nodes 2, 3, and 4. For each of the four nodes, Eqs. (65 1 and (71) are then used to calculate the 
differential cross section for quasiparticles in the vicinity of that node. The coefhcients and differential cross sections 
for kR — 8 are shown in Fig. [3] (Left). Polar plots of the four-node differential cross sections for integer kR from 1 to 
15 are shown in Fig. [3] (Right). In all cases, da /dip has a two peak structure and vanishes for ip — -k. For small kR, 
the peaks are centered about ±7r/2. As kR increases, the magnitude of da /dip increases while the peaks sweep closer 
to the forward direction. 

The single vortex cross section discussed above was calculated to all orders in the linearized potential . However, 
to gain some physical insight regarding the calculated form of our results, it is instructive to consider the more 
transparent Born-limit calculation (valid to first order in the potential) whereby da /dip is proportional to the square 
of the Fourier transform of the scattering potential. For the linearized Hamiltonian, our effective potential is ~ 
—Vf{\/r—r/R^)0{R—r) sin where the 0-function imposes a cutoff at r = i?. Given this input, the Born approximation 
yields an analytical expression for the differential cross section. Details of this calculation is given in Appendix]^ We 
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Figure 3: (Left) Matched coefficients and differential cross sections for kR = 8. The coefficients, b„, denote the contributions 
from eigenstates with total angular momentum n + 1/2. These coefficients and the corresponding differential cross sections 

arc plotted for quasiparticles about nodes 1, 2, 3, and 4. The final row contains the four-node average da/d(p plotted in both 
Cartesian and polar form. The schematic above indicates the directions from which quasiparticles about each node approach 
the non-central effective potential. (Right) Fixed-scale polar plots of the four-node differential cross section for integer kR from 
1 to 15. Note that the cross section magnitude grows with increasing kR. 




Figure 4: Plot of cr|| (the solid line is an interpolation of the numerical data in which the constant large-kR behavior is 
extrapolated to larger kR) 



find 



da 

dip 



IT cos2(v3/2) / _ Ji(2fci?sin(¥^/2)) Y 
Born 8A;1-cos(pV kRsin{ip/2) ) 



(91) 



where Ji (a;) is the first order Bessel function of the first kind. This form comes about because the two-dimensional 
Fourier transform of 1/r yields 1/q where (/^ = |k — k'p = 2A;(1 — cos(^). Hence the square is proportional to 
1/(1 — cos(/)) and goes like l/f"^ for small angles. The rest of the expression results from our cutoff at the vortex edge 
(r — R) and serves to suppress da /dip for ip < 1/kR. Integrating over angles therefore yields a total cross section 
that is linear in kR for large kR. Recall that the transport cross section includes an additional weighting factor of 
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(1 — cos(p). This precisely cancels the leading functional dependence of da /dtp and results in a transport cross section 
that approaches a constant for large kR. Our numerical results for a\\ (see Fig. |4]) are therefore quite logical. Thus, 
the Born approximation is sufhcient to capture the qualitative features of our transport cross sections. 



XII. CONCLUSIONS 



In this paper, we considered the scattering of thermally-excited quasiparticles from a single magnetic vortex in a 
d-wave superconductor. The scattering effect of the vortex can be divided into two contributions, one due to the 
superflow circulating about the vortex and another due to the Berry phase acquired by the quasiparticle upon circling 
the vortex. Both effects are present in the Bogoliubov-de Gennes equation for a single vortex. By applying a 
singular gauge transformation, we have explicitly separated them, the superflow resulting in an effective non-central 
potential and the Berry phase generating antiperiodic boundary conditions for the quasiparticle wave function. In 
a separate papei^S, we have studied the Berry phase contribution in the absence of the superflow effect, essentially 
the Aharonov-Bohm scattering of a massless Dirac fermion. Here, by adopting periodic (rather than antiperiodic) 
boundary conditions, we have considered the superflow contribution in the absence of the Berry phase effect. The 
resulting problem is that of a massless Dirac fermion scattering from a non-central potential (non-central because the 
gap node about which we linearize is shifted away from the origin of momentum space) . We solved this problem by 
calculating the eigenstates inside and outside the vortex, building an incoming plane wave and an outgoing radial 
wave, and computing the resulting differential cross section. We plotted the size and shape evolution of the differential 
cross section as a function of the ratio of the magnetic length R to the de Broglie wavelength 1/k. Intregrating, we 
computed the transport cross section and showed that it saturates to a constant value for large kR. This result 
was previously used in Ref. [13] to compute thermal transport coefficients. Here, we have provided the details of 
its derivation. We have also provided an alternate derivation, and a closed-form solution, valid within the Born 
approximation. These results, taken together with the results of Ref. 12, shed signiflcant light on the nature of the 
scattering of massless Dirac quasiparticles from magnetic vortices. What remains is a complete analysis, taking into 
account both the superflow and Berry phase contributions and the interference between them. 
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Appendix A: Born scattering for massless Dirac equation 



In this appendix, we derive a closed-form expression for the superflow contribution to the single vortex cross section 
within the Born approximatior|22l We begin with a brief recap of Born scattering. The Hamiltonian can be written 
as 

H = H^ + V (Al) 

where Hq is the Kinetic energy operator 

If we denote | > as the energy eigenket of Hq then we have 

H^\4>>=E\cj)> (A3) 

The Schrodinger equation we need to solve is 

{Ho + V)\i,>^E\il}> (A4) 
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The desired solution has the form 

|V>= {E-Ho)-^V\i,> + \^> 
In the position basis one could write it conveniently as 

<x\'ip >=< a: I (/) > + J (fx' <x\iE- HoY^ \ x' ><x' \ V \ x> 

We will express the operator {E — Hq)^^ in the position basis as 



where 



Also note that 



<x\{E- Ho)-^ \x' >= / dVrfV <x\p' ><p' \{E- Ho)-^ I p" >< p" I x' > 



<p' \p >= / <fx' < p' \ X >< x \ p >= 5{p - p') 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



Now we apply this Born scattering method to the single vortex scattering of quasiparticles satisfying the Bogoliubov 
de-Gennes equation 



i7o = Vf{T^Px + Tipy) with E^Vfk 



E ~ Ho — 

(E-Ho)-' = 

<p'\{E- Ho)-' \p"> = 

< X I {E-HoY' \x' > - 1 



Vf{k - r,p, - r,py) = Vf ^ Py^^ 

1 1 f k+p^ Py \ _ ^ (k + T3p.j., + Tipy) 



Vf k"^ -pI-pI \ Py k-p^ 
l (fc + r3P. + np,) ^ 
Vf k"^ ~ p^ 

d^P' . p'-iS-S')i±±j3P^r+riP^ 



Vf 



k^ — p^ 



Vf 



— p'^ 



(AlO) 
(All) 

(A12) 

(A13) 

(A14) 



To evaluate the above integral Eq. (A14l we define the angles specifying the incoming and outgoing momentum 
directions. 



p' = {p',9'), f-f = (I |,a), ip' = a-9' 



(A15) 



Using the above definitions and plugging into Eq. (A14l we have 

1 



<x\iE-Ho) 



f [ , , tp'\x-x'\cosv' (fc + P'jT^ cos(a - + Ti sin(a - ip'))) 
Vf{2nY J P P J k^~p'^ 

-TT 



OO TT 



Vfi2 



I f , f f , f tp'\x-x'\cos v' P'jk + p' COS y>'(T3 cos(q;) + Ti sin(Q;)) 

27r)2 J '^P J '^^^ p_p,2 



dp'-_ 



W/(27r)2 J ^ fc2 -]5'2 




{{kh+p'Mo^h) 
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where Ji = J^^dip'e^"''^°^^ = 27rJo(a) with a = p' \ x — x' \ (J„ being the Bessels functions of first kind) and 
Ma = T3 cos(q;) + Ti sin(Q!). The expression for I2 is given by 



I2 = I Ve'"™"^ cos<p' 



(A16) 



/ Ve^"'°''"' - -*27r£ Jo(a) 



da 

= 2iTiJi{a) 



(A17) 
(A18) 



This enables us to easily rewrite Eq. (AI6I as 



<x\{E- Hoy I x' > 



dp'- 



Vf{2Tr) J ^ fc2 -p'2 



:[{kJo{p' \ x-x' \) + ip'MaJi(p' \ x - x' 



(A19) 



We redefine the variables as u = p'/k , du = dp' /k, and z ^ k \ x — x' \. We also regularize the inverse operator by 
introducing k ^ k + ie with e — !■ 0. The result is 



<a;| {E-Ho + iey^ \ x' > 



k 



du- [Jq{zu) + iu MaJi{zu)] 

1 — + ie 



Vf{27T) 

We remind the reader of the following useful identities: 



du- 



[Jq{zu) +iuMaJi{zu)] 



(A20) 
(A21) 



OC 

/ 



du , , ' = b^K.ib z) for z > 0, Re[b] > , -1< Re[v] < ^ 

+ 0^ 2 

KAz) = i^e'^^'H^z e^t") /or - ^ < Arg{z) < | 



(A22) 
(A23) 



where H^{z ) are the Hankel functions of first kind. Using the above identities in the integrals in Eq. ( A21 1 we obtain 
in the limit e — > 



du 



/ du- 
J u 



uJ^{zu) 
+ (e — i)^ 

v?Ji,{zu) 
2 + (e-i)2 



K,{^i z) = i-Ho{z) 



(A24) 
(A25) 



Hence. Eq. (A21) reduces to 



<x \ {E-Ho + iey^ I x' > 



[Ho{z) + 1 MaHi{z)] 



(A26) 



Now we return to the main idea of setting up Born scattering in the form of self consistent integral equations. 

<x\il)+ >=< a; I > + J d^x' < x \ {E - Hq + ie)-^ \ x' >< x' \ V \ > (A27) 

We apply the eikonal approximation to the above integral equation by assuming | a; |^| x' | which yields | x — x' |~ 
r — f ■ x' ov k' — k f = (k,(j)). In other words a ~ c^, which implies that we are calculating Born scattering to the 
first order only. Now using this approximation, we redefine our variables as Ma — — T3 cos(0) + ti sin(^) and 
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z = k I X — x' \= kr ~ k' ■ x' . Under these conditions Eq. (A26l further simpHfies to, 

<x\{E-Ho + ie)-^ I x' >= [i/o(fer -k' ■x') + i M^Hi{kr - k' ■ x') 



(A28) 



In the Born Hmit, kr ^ 1, we can use the asymptotic form of the Hankel functions of first kind (Hn{z) 
_2_g«(z-f -iif )^ in the above expression to give us 



< X I {E-Ho + ie)-^ I x' 



Ike '2 



[1 + M^y^"- e-'^'-^' 



(A29) 



cos : 



The incoming plane wave in the position basis can be written a& < x \ <f> >= e l ■ ] where k = {k,9) and 
I r 1 = 1 a; |, as in Eq. (53). Using Eq. (A29) we have 

<a;|V+> = <x\(l)> + J (fx'cfx" < x \ {E ~ Hq + ie)-^ \ x' >< x' \ V \ x" >< x" \ (j) > 



k-x I cos 

sin \ 



I k e 2 



'^"^ I d'^x'd'^x" < x' I U I x" > [1 + M^] 



cos : 

sin -, 



The most general form of potential can be written as U = Vq + ViTi + V2T2 + V3T3 = ^nTn where r„ are the Pauli 
matrices and tq = /. In this general representation for V we have 



'/c e * 2 



vji8TT)r 



+ W ^27^6^"" / cPx'cfx"^'^--"-^''-^' J2<x'\Vn\x">[l+ M^]r„ 



COS 

sin 



We now simplify the matrices acting on the spinor 



[1+ Af0]T„ = 



1 + cos 4> sin (j) 
sin d> 1 — cos ( 



1 + cos (j) sin (f) 
sin d) 1 — cos ( 



2/3„(0,0)( 

* sm ^ 



where 



/5o=cos(^— ), /3i =sin(^— ) 
/32 = ism(— — ), /33 =cos(— — ) 



We can write the complete wave function in a suggestive form representing plane wave and scattered wave 

^/r \ sm I 



$(r) 



i k-x 



COS : 

sin -. 



where 



'fee * 2 
1-2 (27r) 



(A30) 
(A31) 

(A32) 
(A33) 

(A34) 
(A35) 



Here f{<j), 9) is the scattering amplitude from which we can easily calculate the differential cross section by 
da 



Ike '2 
vj{2n) 



d^x'd^x"e'^-''"-'^-^' ^<X' \ Vn\x" > Pni(t>,( 



(A36) 
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If V is local and scalar we have the following properties, < x' | Kt | x" >= V{x')S{x' — x") and Vi = V2 = V3 = 0. 
Under these conditions the scattering amplitude simplifies to 



/M = 



6 2 



v^J2tt) 



cos I / d^r e' « 



(A37) 



where (p ^ cj) — 9, q ^ k — k' and q ^ 2k sin (| fc | = | fc' |) and hence 



d(7 
d4> 



I k e ^ 2 ip 

n ,^ ^ COS — 

vj{2n) 2 



6 2 ip ~ 



vV2Tr) 2 



cos -V{q) 



(A38) 



The only part that needs to be evaluated is the integral / d^r e^'^'^V{r) — V {q). The effective potential for the single 
vortex superflow, Eq. (l8| is 



V{r) = VfPsx = -VfPs{r)s'm(j)' 



(A39) 



where Ps{r) — ^(^ — ■^)Q{R — r) and sin0' = sin((^' + 9'). It is important to note the following definitions of the 
angles involved in the scattering. 

k = (fc, 9) (incoming momentum vector) 
k' = (fc, (f)') (scattered momentum vector) 

q = k-k' = {q,9'), f= [r,^') 
If' = <j)' -9', Lp^(j}-9 

Now we are ready to evaluate J d^r e'"^''''V{'r) = V (q) 

V (q) = -Vf [ rdrPsir) [ V e'"''™"'^' sin((p' + 0') (A40) 



= -Vf I rdr^{^~ ^)e{R-r) / V e'"' cos 95' sin 



(A41) 



Using some Bessel functions identities we can reduce this integral to a very compact form involving Bessel function 
Ji. 



V (q) = -invf sine' -{I ~ -^Ji{qR)) 
q qH 

q = 2k \ sm(p/2 \ 



Plugging Eq. (A42l back in Eq. (A38l we get 



da 

dip 



= I ^Vf\l "'Jfo'\ ^^^"^ sin6l'^(l - -^Ji(gi?)) 



1 . 1 



vj{2TT) 



qR 



Also, note the following relations which will help simplifying Eq. (A44l further, 

ip . ip sm9 - sm<j> 

cos — sm9 = cos — ( — ) = cos — ( — — -— - ) 

2 2q' 2 ^2 I sin(/5/2 I ^ 

= —^{cos{(p + 9)+cos{9))sgn{(p) 



(A42) 
(A43) 



(A44) 

(A45) 
(A46) 
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(A47) 



Therefore, simplifying Eq. (A44l further we can write the full cross section as 

da TT 1 , Ji(2fci? I sinG9/2 I) , , ^, 

-r = ^^^7^ 1 hr-r- 7^ '^os O/' + ^ + cos 2 

dip 32fcsin (|) kR | sini^/2 | 

The above quasiparticle scattering cross section is the contribution from one node only. We will now calculate the 
contributions from the other three nodes see (Eq. (90 1). 9 and ip are the Node 1 angles and we calculate the rest of 
the node contribution with respect to this node. 

For Node 2 {02 = -{9 - f ), ip2 = -f) 

oi + 6*1) + cos(6'i) = sin((/7 + 9)+ sin{9) 



C0s((y9l 



For Node 3 (6I3 = (6* + tt), ip'i ^ if) 

cos((/73 + 6*3) + cos(6'i) = — cos((p + 0) — cos(6') 

+ f), ^i = -ip) 

cos((^4 + 9i) + cos((?4) = — sin((y5 + 9) — sin{9) 



For Node 4 {9^ = 



(A48) 
(A49) 
(A50) 



Now we will perform the four node average of the differential scattering cross section, 

da l-^.da. n 1 . Ji{2kR |sinc9/2|),n , , ^, 

- 4g(d^)^-m^(^- kR |sin^/2| '' f<i-<^ + ^)^-os[e)f>-oa.-a.„ 



It turns out that, 



which gives us 



< (C0s(<y9 + 9)+ C0s{9)y > node-average^ 2 COS^ ( - ) (A51) 



da _ cos^(f)^^ Ji(2fcfl |sin(p/2|) ^, ^^^2) 



with our exact results in the weak potential limit. 



dif 16fc sin^d) ' kR I sin(p/2 

'\r 

dip 



We have therefore obtained a closed-form expression for from this straightforward Born-limit calculation. It agrees 
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